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In optical detection setups to measure the deflection of micro-cantilevers, part of the sensing light 
is absorbed, heating the mechanical probe. We present experimental evidences of a frequency shift 
of the resonant modes of a cantilever when the light power of the optical measurement set-up is 
increased. This frequency shift is a signature of the temperature rise, and presents a dependence on 
the mode number. An analytical model is derived to take into account the temperature profile along 
the cantilever, it shows that the frequency shifts are given by an average of the profile weighted by 
the local curvature for each resonant mode. We apply this framework to measurements in vacuum 
and demonstrate that huge temperatures can be reached with moderate light intensities: a thousand 
°C with little more than 10 mW. We finally present some insight into the physical phenomena when 
the cantilever is in air instead of vacuum. 

I. INTRODUCTION 

Silicon micro-cantilevers have today applications in 
many areas of measurement, as chemical and biological 
sensors [1], mass detectors [2, 3], flow meters [4, 5], or 
force sensors [6]. Thanks to their industrial production 
with tight tolerances, and the ability to integrate a sharp 
tip in the fabrication process, they are for example ubiq¬ 
uitous in atomic force microscopy (AFM) [7-10]. The 
detection of their deflection is usually performed opti¬ 
cally: optical lever scheme [11] in most of commercial 
systems, interferometric measurements in the most pre¬ 
cise instruments [6, 12-20]. The former approach is very 
sensitive [21] while remaining simple to set up, and the 
latter reaches the highest resolution in challenging ex¬ 
periments. The limit to the sensitivity of these optical 
systems is due to the shot-noise of the photodetectors: 
this white noise sets the lower limit of detectable deflec¬ 
tions [21, 22]. Its level is linked to the intensity of light on 
the photodetectors: higher intensities translates into bet¬ 
ter signal to noise ratio. Using high power light sources 
is thus appealing. 

The light focused on the micro-cantilever may however 
have some influence on its mechanical response: absorp¬ 
tion of the radiation heats locally the cantilever [2, 23- 
28]. Very high temperatures (several hundred of degrees) 
have been measured for example by Raman scattering 
for laser powers of the order of 10 mW [25, 27, 28]. Per¬ 
manent surface modification [29] or unwanted damage 
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to the tip [23] have been observed in some experiments, 
pointing at the temperature reaching the melting point 
of the material of the cantilever, above 1000°C. The 
effect will be more pronounced in vacuum noteworthily, 
since heat has no other efficient way out than through the 
tiny cross-section of the cantilever, producing a notice¬ 
able temperature rise of the probe. This photo-thermal 
effect can be useful in some applications, for example to 
drive a cantilever at its resonance in dynamic measure¬ 
ments [23, 24, 30-32], or to heat the AFM tip and control 
local phase changes of a sample [33]. On the contrary, 
this can be an issue for temperature sensitive samples, or 
for low temperature experiments [6, 18], those aiming at 
reducing thermal noise for example. 

Silicon is a poor mirror for visible light, with only 37 % 
of reflectivity, thus radiation absorption is large in raw sil¬ 
icon micro-cantilevers [34]. Reflex coatings (Aluminium, 
Platinum, Gold) are often used to minimize this issue, 
but they add some viscoelastic damping to the probe me¬ 
chanical response, lowering its quality factor and thus its 
performance in dynamic operation [35-38]. Understand¬ 
ing the effect of heating a cantilever through the mea¬ 
surement system is thus important for its operation [39]. 
In this article, we present an approach to quantify the 
temperature rise and its consequences on the mechanical 
response of micro-cantilevers, in vacuum and in air. The 
analysis applies consistently across all flexural modes of 
the cantilever, a welcomed point in nowadays growing 
scope of application of multifrequency AFM [10]. 

The article is organized as follows: in part II, we 
present experimental evidence of a frequency shift of the 
resonant modes of a cantilever when the light power of 
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the optical measurement set-up is increased. This fre¬ 
quency shift is a signature of the temperature rise [26, 27, 
39, 40], but its dependence on the mode number demon¬ 
strates that this temperature cannot be uniform on the 
cantilever. In part III, we present an analytical model to 
take into account the temperature profile along the can¬ 
tilever, and show that the frequency shifts are given by an 
average of the profile weighted by the local curvature for 
each resonant mode. In part IV, we apply this framework 
to measurements in vacuum. We demonstrate that huge 
temperatures can be reached with moderate light inten¬ 
sities: a thousand °C with little more than 10 mW. In 
part V, we finally present some insight into the physical 
phenomena at play when the cantilever is in air instead of 
vacuum, before giving a general conclusion to this work. 


II. EXPERIMENTAL DATA: FREQUENCY 
SHIFT 


We perform the experiments on two different geome¬ 
tries of cantilevers: cantilever ClOO is L = 500 pm long, 
W = 95 pm wide, and ff = 0.75 pm thick (NanoWorld 
Arrow TL8), while cantilever C30 is L = 500pm long, 
W = 30 pm wide, and H = 2.8 pm thick (BudgetSensors 
AIO-TL). Both are uncoated tipless AFM silicon can¬ 
tilevers. Geometrical dimensions were measured using 
a scanning electron microscope (SEM), with uncertain¬ 
ties around 1% for L and W, and 5% for H. The ther¬ 
mal noise measurements are performed using a quadra¬ 
ture phase differential interferometer [20], with the sens¬ 
ing laser beam focused on the free end of the cantilever, 
and the reference beam on the silicon chip where the 
cantilever is clamped (see inset of figure I). We use a 
stabilized solid state laser from Spectra Physics, with a 
100 mW maximum output power at 532 nm. We couple 
this laser through an optical fiber into our measurement 
set-up, and an intensity / up to 25 mW can be used for 
the sensing beam. This intensity can be tuned to lower 
levels using neutral density filters at the laser output, en¬ 
suring highly stable light intensity and frequency during 
the whole measurement process. 

We record the power spectral density (PSD) of the 
brownian deflection of the cantilever using our interfer¬ 
ometer, and plot the results in figure 1. Such a spec¬ 
trum is usually computed by an average on 20 PSD, 
each evaluated on a 5 s deflection measurement sampled 
at 1.4 MHz. The 6 first normal modes of the cantilever 
ClOO are clearly seen of this spectra. A few higher or¬ 
der modes could also be measured, but we limited our 
study to these 6 modes. To ensure precise positioning of 
the sensing laser beam on the cantilever and avoid any 
spurious effect from thermal drifts, we actually use the 
7th resonant mode of this cantilever : we look for the last 
spatial node of this mode close to the free end of the can¬ 
tilever. Killing any thermal noise in the spectrum at the 
corresponding resonance frequency defines the longitudi¬ 
nal position of the sensing beam with better than 3 pm 
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FIG. 1. (Color online) — Thermal noise measurement on 
cantilever ClOO: we plot the power spectrum density (PSD) 
of thermal noise induced deflection as a function of frequency, 
for various incident light power I on the lever (I = 0.7 mW 
to ll.OmW), in vacuum (P = 4Pa). The flat base line of the 
measurement is due to the shot noise on the photodetectors, 
it decreases when I increases. All the resonance frequencies of 
the cantilever decrease when I grows, as illustrated in the two 
bottom figures for modes 1 and 3. The inset of the top figure 
illustrates the principle of the measurement: the fluctuations 
of deflection are recorded through the interference of two laser 
beams, one reflected on the cantilever free end, the other on 
the chip holding the cantilever [20]. 


repeatability. Centering the laser spot on the cantilever 
is also done with the same precision killing any thermal 
noise corresponding to torsional modes of the lever. A 
similar procedure can be applied to cantilever C30, for 
which we characterize up to 4 modes. 

The background noise of the instrument is as low as 
5 X 10“^®m^/Hz, as illustrated by the noise floor at the 
highest light intensity measurement in figure I. The main 
contribution to this background noise is the shot noise of 
the photodiodes used in our optical setup, it thus de¬ 
creases when light intensity increases. It is therefore 
tempting to use the highest light power available to in¬ 
crease the dynamic of our measurement. We were able 
to lower the background noise down to 10“^° m^/Hz for 
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some cantilevers with a reflective coating. However, the 
sensing light is also partially absorbed by the cantilever, 
especially for raw silicon cantilevers presented in this arti¬ 
cle, leading to a heating of the sensor, potentially modify¬ 
ing its response or integrity. And indeed, a closer look on 
the resonant peaks of the PSDs in figure 1 demonstrate a 
light power dependence: the resonance frequencies shift 
when / changes [26, 27, 39]. To study this effect, we fit 
each of these resonances with a lorentzian curve to ex¬ 
tract the resonance frequency /„ for each mode number 


We plot in figures 2 and 3 the measured relative fre¬ 
quency shift Sfn/fn for each mode number n as a func¬ 
tion of the light intensity I of the sensing laser beam at 
the cantilever end. /° is defined as the ordinate at the ori¬ 
gin of a linear fit of /„ versus / for low laser powers, such 
that Sfn = 0 for / = 0. These curves depend on mode 
number, on the pressure of air around the cantilever, and 
on the cantilever shape. We also observe some nonlinear 
behaviors of the frequency shift as a function of light in¬ 
tensities for the larger cantilever ClOO when / > 6 mW. 
The laser spot size is tuned around 10 pm in diameter 
in these measurements to avoid most of light spilling out 
of the cantilever. The non-linear effects tend to occur 
sooner for smaller spot sizes. 

A simple argument can be used to understand the 
global behavior of the resonance frequency [26, 27, 39, 
40]: each mode can be pictured as a harmonic oscillator 
of mass TOeff, whose stiffness kn is proportional to the 
Young’s modulus E of silicon (the cantilever material). 
The increase of temperature, due to light absorption, in¬ 
duces a softening of the cantilever: the temperature co¬ 
efficient of E is negative [39, 41-45] 
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The resonance frequencies of the normal modes should 
therefore decrease as light intensity increases: 
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However, this simple argument does not explain why the 
frequency shift is mode number dependent (and even less 
how the shift can be positive in air for cantilever ClOO). 
In the next part, we show that the mode number de¬ 
pendence can be explained by the coupling between the 
spatial mode shape of the normal modes and the temper¬ 
ature profile in the cantilever. 


III. ANALYTICAL DESCRIPTION: 
FREQUENCY SHIFT OF RESONANT MODES 

In the Euler-Bernoulli framework to describe the 
micrometer-sized mechanical beam, we assume that the 


FIG. 2. (Color online) — Cantilever ClOO: relative frequency 
shift Sf„/fn for modes n = 1 to 6 as a function of the in¬ 
cident light power I on the lever, at atmospheric pressure 
P — 10® Pa, top) and in vacuum (P = 65 Pa middle, P = 4 Pa 
bottom). The dashed lines are linear fits on the initial part 
of the curves (7 < 6mW). The curves shapes are nonlinear, 
and depends strongly on pressure and mode number. 


cantilever length L is much larger than its width W, 
which itself is much larger than its thickness 77. The 
flexural modes of the cantilever are supposed to be only 
perpendicular to its length and uniform across its width. 
The deformations can thus be described by the deflec¬ 
tion d{x, 7), with t the time, and x the spatial coordinate 
along the beam. The equation of motion for the can- 
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FIG. 3. (Color online) — Cantilever C30: relative frequency 
shift 6f„/fn versus incident light power I on the lever, at 
atmospheric pressure P = 1 x 10® Pa, top) and in vacuum 
(P = 10Pa bottom). The dashed lines are linear fits. 


tilever reads [46, 47]: 




B(x)/(«)g 




( 4 ) 

where /r is the mass per unit length, El the bending 
stiffness {E Young’s modulus, I the second moment of 
inertia) of the cantilever, t) the force per unit 

length due to the surrounding fluid (if any) acting on the 
cantilever, and /®’'*(a;,t) the external force. 

We model the influence of the fluid according to Sader’s 
model [47]: t) is supposed to be linear in the local 

deflection d{x,t). We therefore write in Fourier space: 


/^“'^(a;, uj) = ^fUJ^r(uj, iy)d(x, u) (5) 


where fJ.f = 7rpVF^/4 is the mass per unit length of a 
cylinder of diameter W of fluid of density p, v is the 
fluid kinematic viscosity, and F is the complex hydro- 
dynamic function of a rectangular cantilever [47]. The 
real part F,. of F translates into inertial effects, while the 
imaginary part F^ accounts for the viscous damping. In 
vacuum, dissipation occurs only through internal mecha¬ 
nisms (viscoelastic damping) in the cantilever [36, 48]. It 


can be modeled using a complex Young’s modulus, but 
for simplicity we shall mainly consider viscous dissipation 
in this paper. 

Due to the heating of the cantilever by the sensing 
laser beam, its temperature is not uniform. In line with 
previous geometrical hypotheses, we will suppose that 
the temperature profile is uniform across the section of 
the cantilever and time independent, thus of the form 


T{x) = To + e{x) (6) 


where Tq is the clamping base temperature (we make the 
hypothesis that the clamping base is at room temper¬ 
ature thanks to its macroscopic size. d{x) is thus the 
temperature difference profile along the cantilever. This 
temperature profile will locally change the physical prop¬ 
erties of the cantilever and surrounding atmosphere if not 
in vacuum. As a first approximation for the cantilever, we 
will only consider the variation of Young’s modulus and 
dissipation with temperature, and neglect any change in 
the physical dimensions. Indeed, the coefficient of lin¬ 
ear thermal expansion of silicon ai ~ 2.6 x I0“®K“^ is 
one order of magnitude smaller than the temperature co¬ 
efficient of its Young’s modulus as ~ —64 x I0“®K“^ 
(numerical values at room temperature [42]). Moreover, 
we will only consider first order correction for the Young’s 
modulus: E{T) = To(I -I- a^^). Note that we only need 
to assume that <C I (which is a light assumption 
since it means 0 <S I/as = 15 000 K), but not that 
0 <C Tq. We therefore rewrite equation (4) in Fourier 
space in time to explicit the temperature-dependence of 
the bending stiffness: 


fcp d'^ 
3 dx^ 


, ..d^d 
{l + aEe{x))^ 


— {m -I- m/F(w, v))ui^d = Lp^^{x,uj) 


( 7 ) 


where m and fep ^re the cantilever mass and static stiff¬ 
ness (fcp is defined at room temperature Tq), mf = Lfif, 
and X is now normalized to the length L {x = 1 at the 
free end of the cantilever). In absence of external forc¬ 
ing (normal modes), the boundary conditions for Eq. (7) 
correspond to the classical clamped {x = 0) and free end 
(a: = I) conditions: 


d{x=0, oj) = 0 


d'^d{x,uj) 


dx^ 


= 0 


21 = 1 


dd{x, uj) 
dx 

d^d{x,oj) 

dx^ 


x—0 

a;=l 


= 0 

= 0 


( 8 ) 


A. Naive approach : uniform temperature 

Let us first study the case where 9{x) = 0. We in¬ 
troduce the operator Cq = d'^jdx'^ in the space V of 
functions matching boundary conditions (8). This oper¬ 
ator is self-adjoint on this space for the integral scalar 
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product : 


(/) © '0 



(f>{x)'tp{x)dx —>■ £o0 0 i> 


4> 0 Coip 


(9) 


The eigenvectors of operator Cq form an orthonormal 
basis of V. They are the classic normal modes of the 
Euler-Bernouilli model of the free-clamped mechanical 
beam: 


0° (x) = (cosa„x — cosha^a::) + i3„(sina„a: — sinha^x) 
with 

^ cos an + cosh a„ 

" sin Qfn + sinh an 

and where a„ is the solution of equation 

1 -f cos an cosh a„ = 0 

which leads to ai = 1.875, 02 = 4.694, ..., and Q!„ = 
(n — l/2)7r for large n. The eigenvalues A° of the operator 
£0 corresponding to each normal mode are simply A° = 
Otn- ^0(l^n = K.4‘n = “n'/'n- 

We project equation (7) on this basis to get the evolu¬ 
tion of the amplitude of each mode: 


( 10 ) 

( 11 ) 

( 12 ) 


2 I ^0 4 I • 


d„(w) = £r‘(w) (13) 


with dn = dQ(j)n^ = L/®’'*00°, 7eff = wm/ri(a;,i^) 

and Weff = m + ■mfTr{uj,v). T(u},v) is usually slowly 
varying with w, and dissipation is low (large quality fac¬ 
tor) , thus for every resonance a simple viscous damping is 
a good approximation. The conclusion would be the same 
with a viscoelastic damping, though the damping coeffi¬ 
cient 7 eff would present an explicit mode dependence. We 
therefore have a collection of independent quasi-harmonic 
oscillators, of effective mass meff, stiffness = a^fco/S, 
and damping coefficient 7 eff. 

The resonant pulsations of the oscillators are linked 
to the spatial eigenvalues by the dispersion relation: 

2 ^0 4 ^0,0 /T A \ 

= ya„ = yA„ (14) 

If we now consider a uniform temperature rise 6 = AT, 
assuming meff to be constant, we only need to change the 
value of fco by fco(l + aEAr). The relative frequency shift 
is thus independent of the mode number n and evaluates 
to 


SuJn l(5fco 1 .rr 

(-^n ^ /Cq ^ 


(15) 


cantilever : asdix) <C 1. We introduce the operator 
£ = 5^/5x^[(l -I- aE9{x))d‘^/dx^], which is still self- 
adjoint on space T>. The eigenvectors 0„ associated to 
the eigenvalues A„ of this operator form an orthonormal 
basis of V. They are defined by: 

d^4>n = An0n (16) 


Let us expand this relation to first order in ue with £ = 
£o -I- q;_e£i, 0 „ = 0 ° -I- a_E 0 i( and A„ = A° -I- aE^h- 

£o0i+£i0° =A°0i+Ai0° (17) 

Projecting this last equation on 0° and using the prop¬ 
erty of £o being self-adjoint, we find 

Ai=0O0£i0O (18) 

Using the explicit expression £i = jdx'^[9{x)d'^ Jdx^] 
and some simple calculations, it is straightforward to 
show that 

= (19) 

For any given temperature prohle 9(x), it is therefore 

easy to compute the induced shift of the spatial eigen¬ 
value. We can then use the dispersion relation (14) to 
deduce the shift in resonance frequency: 


SijJn 

iUji 


1 aE-^n 

2^ 

1 1 

/ dx9(x) 


\ dx^ J 


( 20 ) 

( 21 ) 


The frequency shift is thus mode dependent: the soften¬ 
ing (as < 0) of the cantilever due to the temperature 
rise is weighted by the square of the local curvature, that 
is where the cantilever is most bended. Note that if we 
choose a constant temperature profile 9 = AT, we re¬ 
cover the mode independent frequency shift of equation 

(15).^ 

It is important to notice that the formal expression 
(21) of the frequency shift is independent of the dissi¬ 
pation mechanism, whether viscous (from the surround¬ 
ing fluid) or viscoelastic (internal). The presence of air 
around the cantilever only influences the frequency shift 
due to inertial effects (variation of meff with T), and 
through its thermal conductivity which alters the tem¬ 
perature profile. 


C. Application: linear temperature profile 


B. Arbitrary temperature profile 

We now allow any temperature profile, but work in 
the limit of small perturbations to the stiffness of the 


Let us consider the case of a linear temperature pro¬ 
file 9{x) = xAT, where AT is the temperature increase 
at the free end of the cantilever. We report in table I 
the expected coefficients linking the relative frequency 
shift of this temperature profile to the relative frequency 
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shift of a cantilever with uniform temperature AT. This 
coefficient is close to 1/2 for large mode numbers which 
present a strong curvature all cantilever long: we observe 
the effect of the mean temperature ATj2. For low mode 
numbers however (especially mode 1), this coefficient is 
smaller than 1/2: curvature is larger close to the clamped 
end which is at lower temperature than average, thus the 
frequency shift is not that large. 


6mW for cantilever ClOO to avoid non linearities). The 
slope of these fits, <5/„//„/, are plotted against in fig¬ 
ure 4. For both cantilevers in vacuum, the mode depen¬ 
dence is correctly captured by our analysis. The prefactor 
is very similar for both cantilevers: though their geom¬ 
etry is different, they share very similar length L, cross 
section WH, and coefficient of absorption of light by the 
cantilever a. 


Mode number 

1 

2 

3 

4 

5 

6 

dujji 2i 

_QbAT 

0.193 

0.406 

0.476 

0.483 

0.490 

0.493 


TABLE I. Relative frequency shift for a cantilever with a lin¬ 
ear temperature profile normalized to the relative frequency 
shift of a cantilever with the same maximum temperature but 
a flat profile. 


IV. TEMPERATURE PROFILE OF A 
CANTILEVER IN VACUUM HEATED AT ITS 
EXTREMITY 


Let us first consider the case of the measurements in 
vacuum. In this case, the main dissipation process for the 
heat due to the absorption of light by silicon is thermal 
conduction along the cantilever. The stationary temper¬ 
ature profile can be estimated using Fourier law for heat 
transfer: 


J = 


al 

Wh 


\{T{x)) 


1 dr 

L da; 


( 22 ) 


where J is the absolute value of the heat flux (in W/m^), 
a is the fraction of light absorbed by the cantilever, and A 
the thermal conductivity of silicon. We assume here that 
the heat flux is homogeneous across the cantilever cross 
section WH. The temperature profile can be computed 
from this equation. 


A. Low intensities: linear approximation 

For low light intensities /, we neglect the temperature 
dependence of A, the temperature profile is then trivially 
linear: 


e{x) 


Lai 

WHXo 


(23) 


In such a case, the dependence of the frequency shift on 
mode number n should be modeled by the coefficients 
of table I, computed for a linear temperature profile : 


uJn ”2 ^ ^ 


To test this dependence, we perform linear fits of the 
curves of figures 2 and 3 (considering only data below 
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FIG. 4. (Color online) — Slope of the relative frequency shift 
versus light power as a function of coefficient (expected 
mode dependence for a linear temperature profile). The pro¬ 
portionality (dashed lines) is reasonable for both cantilevers 
in vacuum. At atmospheric pressure, a global vertical offset 
(mode independent contribution) can be observed. 


Using Eq. (24), we can also estimate the maximum 
temperature reached at the cantilever free end from the 
experimental frequency shift. As seen in figure 5, this 
estimation is coherent for the different modes. The tem¬ 
perature rise is significant: around 100 K for a 3 mW il¬ 
lumination, and up to BOOK for 12mW ! For the largest 
light intensities, the curve is not linear anymore: the de¬ 
pendence of A on temperature cannot be neglected. 


B. High intensities: non-linear approach 

We plot in the inset of figure 6 the behavior of X(T) 
from references [49, 50] (experimental data and model). 
We can use these data to compute the temperature profile 
as follows: defining 

A(r) = [ A(T')dT', 

4 To 


( 25 ) 
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FIG. 5. (Color online) — Temperature increase at the free 
end of the cantilever estimated from the frequency shift in 
vacuum, assuming a linear temperature profile. This approx¬ 
imation tends to be false (underestimated) at higher power 
where nonlinearities are higher, and the estimation of AT for 
the different modes tends to be more disperse. We see how¬ 
ever that large temperature gradients are already present for 
a few mW illumination. 
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FIG. 6. (Color online) — Theoretical temperature profiles 
computed for a cantilever with a cross section WH = 71 pm^ 
and an absorption of light a = 0.6, taking into account the de¬ 
pendence of thermal conductivity of silicon A on temperature. 
The only thermal transfer process taken into account here is 
conduction along the cantilever. The profile is clearly nonlin¬ 
ear for light intensities above 6mW. Inset: dependence of the 
thermal conductivity of silicon A on temperature T [49, 50]. 
The variation of A is significant and should be taken into ac¬ 
count. 


we have by integrating equation (22) over x that JLx = 
A(T(x)), from which the temperature profile is given by 

0{x) = A-\jLx)-To. (26) 

From the data in the inset of figure 6, we can esti¬ 
mate numerically A(T) and thus what the temperature 
profiles are for various light intensities I. We plot these 
estimations in figure 6 for cantilever CIOO, assuming an 
absorption coefficient a = 0.6. A more accurate esti¬ 
mate of the temperature profile would require to know 
the temperature dependence of a; here we consider for 
simplicity a typical, temperature-independent value of 
a. The nonlinear shape of the temperature profiles are 
clearly visible for light powers above 6mW. We report 
the maximum temperature expected from these non lin¬ 
ear profiles (at the free end of the cantilever) in figure 
7. The agreement with the experimental data is good, 
especially as the measurement points are relying of a lin¬ 
ear temperature profile analysis. Note that in principle, 
non-linear temperature profiles should also be considered 
in the analysis of the frequency shifts. However, since 
the precise, temperature-dependent value of the light ab¬ 
sorption coefficient is unknown, the accuracy of the de¬ 
termination method would be difficult to assess, and we 
thus stick to the analysis based on a linear temperature 
profile. 

Since the melting point of silicon is = 1410 °C, we 


see that this analysis predicts dramatic consequences for 
the cantilever for light intensities higher than 15mW ! 
And indeed, for the highest intensities we have probed 
in vacuum, some cantilever have been molten at their 
extremity. We present for example in figure 7 a SEM 
image of a C30 cantilever, that could be shorten from 
its initial length L = 500 pm to only 275 pm using I = 
25 mW. 

Using the temperature profiles of figure 6, we can esti¬ 
mate the frequency shifts from equation (21). The results 
are plotted in figure 8. The curves are qualitatively very 
similar to the experimental measurements of figure 2. A 
quantitative agreement cannot be achieved anyway, and 
we can point a least 3 levels of approximation : 

1. the temperature dependence of the Young’s modu¬ 
lus of silicon would have to be better known in the 
temperature range of interest. 

2. the coefficient of absorption a depends on temper¬ 
ature. 

3. the temperature is not uniform across the cantilever 
cross section close to the free end where the laser 
is focused. 

The general framework presented here is thus appro¬ 
priate to describe the behavior of cantilevers in vac¬ 
uum heated by light absorption at their free extrem- 
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FIG. 7. (Color online) — Maximum temperature reached at 
the end of cantilever ClOO as a function of light intensity: 
model (plain line) and experimental data from figure 5 (mean 
temperature extracted from the 6 measured modes, the uncer¬ 
tainty corresponds to the standard deviation between those 
6 points). Inset : SEM Image of a C30 cantilever heated by 
a laser power of 7 = 25 mW: the temperature raised above 
= 1410 °C, leading to the melting of silicon. The initial 
length L = 500 pm has been shorten to 275 pm when moving 
the laser focus towards the base. 



FIG. 8. (Color online) — Theoretical (plain lines) and exper¬ 
imental (symbols) relative frequency shift Sfn/fn for mode 
n = 1 to 6 as a function of incident light power I on the 
lever. The theoretical curves have been computed for a can¬ 
tilever with a cross section WH = 71 pm^ in vacuum, with 
an absorption of light of a = 0.6, taking into account the de¬ 
pendence of thermal conductivity of silicon A on temperature. 
The experimental data is that of figure 2 in vacuum at 4 Pa. 


ity. Through the proposed analysis, the frequency shifts 
prove to be a robust way of estimating the temperature. 


V. INERTIAL EFFECTS OF THE 
SURROUNDING FLUID IN AIR 

In air, the observed frequency shift depends strongly 
on the cantilever shape, and can even be opposite to 
the expected trend: the resonance frequency increases 
with laser power for cantilever ClOO. However, addi¬ 
tional effects have to be taken into account due to the 
presence of the atmosphere: the temperature rise of the 
cantilever heats the surrounding air, increasing its viscos¬ 
ity and lowering its density. The latter effect leads to a 
decrease of the inertial effects due to the air moving along 
with the cantilever, hence help raising the resonance fre¬ 
quency [40]. This effect is more pronounced on larger 
cantilevers since the volume of air implied is higher. 

We did not conduct a systematic study of the interac¬ 
tion between this effect and a non uniform temperature 
profile, but a first estimation of the effect of the temper¬ 
ature rise in air can be done with the Sader model [47]. 
Indeed, from equation 14 for a uniform temperature Tq, 
we can compute for every cantilever the resonance fre¬ 
quency shift: 


5uJn _ ISkp 1 Snies ^ 27 ^ 

UJ-fi 2 /cq 2 TTlefi 

The first term accounts for the temperature dependence 
of the Young’s modulus of Silicon, while the second term 
accounts for the variation of the effective mass of the 
oscillator. Let us explicit nies: 


meff(To) = m + ^LW^p{To)rr{uj, iy{To)) (28) 


In this equation, we neglect again the temperature de¬ 
pendence of the geometry {L and W supposed constant), 
but consider explicitly the variation of the air density p 
and kinematic viscosity v with Tq [51]: 


p{Tq) = 

7o 


v{To) = Vi, 


To + S 



(29) 

(30) 


with Pi, = 1.29kg/m3, T* = 273 K, S = 110AK and 
Vi, = 1.32 X 10“^m^/s. From these equations and the 
explicit expression of the hydrodynamic function r(a;, i^) 
in reference [47], we can compute Srues /meg as a function 
of a uniform temperature Tq. 

We plot in figure 9 the frequency shift expected for a 
silicon cantilever as a function of Tq. We observe that 
when Tq increases, the inertial effects due to the decreas¬ 
ing density of air are opposite to the softening of the 
Young’s modulus of silicon. For the wider cantilever, the 
former effect is dominant, whereas for the narrower one, 
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FIG. 9. (Color online) — Theoretical relative frequency shift 
Sfn/fn versus uniform temperature T of cantilever and air 
according to the Sader model. The shift is the sum of two 
opposite effects when temperature rises: the softening of the 
cantilever (only effect in vacuum), and the decreasing of the 
inertial effects due to the surrounding atmosphere. The latter 
effect is dominant for the wider cantilever at low T. For a 
uniform temperature, the mode dependency is weak. 


it is just a correction to the latter. Though the temper¬ 
ature profile is not considered in this model, a similarity 
with the experimental data of figure 2 can be seen : the 
frequency shifts of intermediate modes present a max¬ 
imum and then decay, when inertial effects due to the 
fluid vanish at high temperature. 

If the temperature profile was to be considered, dissi¬ 
pation of heat through air would lead to a temperature 
rise located at the end of the cantilever only. The in¬ 
ertial effects would thus be more important on mode 1, 
for which the highest deflection are close to the end. At 
the same time, the base of the cantilever, where most of 
the deformation occurs, will not see the temperature rise, 
thus softening of the cantilever is negligible. And indeed, 
the maximum for mode 1 is not seen in the experimental 
data of figure 2, the inertial effects being always domi¬ 
nant. 

As a hnal note, if we focus on the frequency shift due 
to the inertial effects only in this model, we checked that 


it depends only weakly on the mode number. The slope 
of the frequency shift versus light power should therefore 
just be translated by a mode independent contribution 
in air with respect to vacuum, as observed on the exper¬ 
imental data of figure 4. 

VI. CONCLUSION 

To measure the temperature profile of the cantilever, 
we propose to use the frequency shifts triggered by the 
softening of silicon upon heating. The proposed model 
adequately describe the experimental data in vacuum be¬ 
low 200 °C, and could be quantitative at higher temper¬ 
ature provided the Young’s modulus of silicon is better 
characterized in the adequate temperature range. Our 
approach is robust for all flexural modes of the cantilever, 
and relates the frequency shift to the average of the tem¬ 
perature profile weighted by the local curvature. Coef¬ 
ficients to apply to modes 1 to 6 are given in table I 
for direct application to linear temperature profiles. As 
for observations in air, a qualitative description of the 
measured frequency shift has been proposed through the 
dependence of inertial effects due to the surrounding fluid 
in its temperature. 

As demonstrated in this article, strong heating of a raw 
silicon AFM cantilever in vacuum can be achieved using 
light absorption of a laser focused at the free end of the 
cantilever. The temperature rise depends on the length 
and cross section, but can be larger than a thousand °C 
for a 500 pm long cantilever with little more than 10 mW 
at 532 nm. Such large laser intensities are not common 
in standard AFM operation, but the usual 1 mW is al¬ 
ready enough to produce a 30K increase, which should 
be considered for temperature sensitive samples. Reflex 
coated cantilevers, or shorter ones, will obviously reduce 
this issue. On the other hand, if one wants to heat locally 
a sample (to trigger phase changes or modify other prop¬ 
erties) , high temperature of the AFM tip can be reached 
through a simple external laser focused on the cantilever 
free end. In this case, the current framework provides an 
easy procedure to estimate the temperature of the AFM 
probe. 
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